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We study a six-dimensional braneworld model with infinite warped extra dimensions in the case 
where the four- dimensional brane is described by a topological vortex of a U(l) symmetry-breaking 
Abelian Higgs model in presence of a negative cosmological constant. A detailed analysis of the 
microscopic parameters leading to a finite volume space-time in the extra dimensions is numerically 
performed. As previously shown, we find that a fine-tuning is required to avoid any kind of singular- 
ity on the brane. We then discuss the stability of the vortex by investigating the scalar part of the 
gauge-invariant perturbations around this fine-tuned configuration. It is found that the hyperstring 
forming Higgs and gauge fields, as well as the background metric warp factors, cannot be perturbed 
at all, whereas transverse modes can be considered stable. The warped space-time structure that 
is imposed around the vortex thus appears severely constrained and cannot generically support 
nonempty universe models. The genericness of our conclusions is discussed; this will shed some light 
on the possibility of describing our space-time as a general six-dimensional warped braneworld. 

PACS numbers: 04.50.-(-h, ll.lO.Kk, 98.80.Cq 



I. INTRODUCTION 

Following the advent of string theory ^ and its im- 
plication that space may have more than the usual three 
dimensions (in the Kaluza-Klein way) came the sugges- 
tion that the extra dimensions could be much larger than 
previously expected, would it be because of a smaller 
value of the Planck energy in the bulk Q, H, Q , or be- 
cause of a large curvature in the (infinite) extra dimen- 
sions H, 1^' A novel idea came into play with the as- 
sumption that we live on a hypersurface, a three-spatial 
dimensional "brane" , embedded in a larger dimensional 
warped space-time bulk 0, |1| . 

For any higher dimensional Universe model, it is es- 
sential to confine gravity since gravitation is experi- 
mentally tested to be three dimensional on many dif- 
ferent scales, ranging from the millimeter to a few 
Mega-parsecs jl^ For a five-dimensional anti-de 

Sitter bulk, gravity was shown to be localized on the 
brane JMHand to lead to a viable cosmological frame- 
work [Hill El El El m El El E provided the 
brane and bulk cosmological constants are adjusted by 
hand. The situation is not yet settled concerning the 
cosmological perturbations induced in the brane |2l| . 
although there are some indications that such models 
should satisfy more stringent constraints than previously 
expected 22] . In the case the brane is modeled as a do- 
main wall-like topological defect, this fine-tuning trans- 
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forms into a tuning of the underlying parameters (masses, 
coupling constant and bulk cosmological constant) 23j. 

Many mechanisms have been proposed to confine the 
other known interactions and their associated particles: 
scalar J4J, gauge bosons 25^ 26 ■, 21, 28, 29, 30, 31, 3i, 
[3^ . and fermions In the latter case the mech- 

anism relies on a generalization of the cosmic string 
case 0, 0, ^ . The fermions are trapped in the brane 
in the form of massless zero modes and some can even 
become massive, although their mass spectrum is not 
compatibl e [221 with the observed one [Sg- It was also 
suggested 39] that the electroweak Higgs field, and thus 
the origin of electroweak symmetry breaking, could be 
understood from the existence of an extra dimension in 
the form of a transverse gauge field component. Most 
of these works were based on the simplifying assumption 
of a reflection symmetry with respect to the brane, al- 
though a more refined treatment, not assuming such a 
symmetry, appears possible IH |4l|, |43| • 

Most of the relevant discussions on braneworld mod- 
els have been restricted to the case of one spatial ex- 
tra dimension, as advocated e.g. in the framework of the 
eleven dimensional realization of theory proposed by 
Hofava and Witten fi^. Moreover, the underlying brane 
model of the Universe is often assumed to be infinitely 
thin in the transverse direction, so that (i) the induced 
gravity stems essentially from Darmois - Israel junction 
conditions [HElEllllElElElElllSand (n) is 
mostly independent of the microstructure of the brane, if 
any. No such general condition is available in the less re- 
strictive situation of more than one extra dimension that 
is the subject of the present work. 

To study braneworlds with more than one extra dimen- 



2 



sion, it is necessary to specify the microstructure of the 
brane to fix a model, taking into account in particular the 
possibly finite thickness of the brane 0|- In particular, 
one needs to regularize long range self-interaction forces 
(including gravity). In the case the force derives from a 
potential (e.g. in linearized gravity), the self- interaction 
potential is well behaved outside the brane but will be 
singular on the brane in the thin brane limit. One needs 
to introduce a UV cutoff associated with the underlying 
microstructure The only case where such a proce- 
dure can be avoided is that of hypermembranes that can 
be satisfactorily treated without recourse to regulariza- 
tion. 

Much work has been devoted to warped geometries in 
six dimensions IM El El El ES IH 113 ■ When con- 
sidering explicit realizations in terms of an underl ying 
topological defect forming field model [H 0, IH [sj 
it seems that the six-dimensional case represents a lim- 
iting situation: for more than two extra dimensions, it 
was found that it is not possible to confine gravity on 
a strictly local brane although global topological defect 
configurations which allow gravity confinement might ex- 
ist |53,|53|- At least two questions arise: are the proper- 
ties of gravity dependent on the microscopic structure of 
the brane and is the chosen microscopic model consistent 
with M— theory? The second question has started to be 
addressed in Ref. and it turns out that, for many 
purposes, it is possible to consider defectlike realizations 
of branes, as in the present article. 

Before considering the possibility of trapping [60l IbH 
particle fields in a hyperstring embedded in an anti-de 
Sitter six-dimensional bulk space-time (adSg ) , it is neces- 
sary to determine the background structure itself within 
a given field content, and decide whether it is possible 
to localize gravity in the Universe thus obtained, thereby 
generalizing the five-dimensional case. This article is de- 
voted to this task and accordingly models the brane by a 
vortex of a U(l) Abelian gauge-Higgs model. By means 
of a numerical exploration of the parameter space, we 
discuss different classes of solution exhibiting an anti-de 
Sitter space-time at infinity. It appears that they are 
generically associated with a conical or curvature singu- 
larity at the brane location, except for some fine-tuning 
between the model parameters. In the regular case, our 
approach agrees with the numerical results of Ref. |62| . 

Assuming this fine-tuning, we then go on to analyze 
the stability of the regular solution by studying the scalar 
modes of the gauge invariant perturbations, as originally 
suggested in Ref. "b^]. Restricting our attention to the 
lowest angular momentum modes, we show that the only 
acceptable perturbations of the nonvanishing background 
quantities, i.e. the hyperstring forming fields and metric 
warp factors, are the vanishing ones. The conditions im- 
posed on these perturbations to be acceptable being to be 
initially finite with respect to the background value away 
from the brane (a condition necessary in order to ensure 
that the bulk is close to anti-de Sitter space), and to be 
well-behaved on the brane. We also study the perturba- 



tions of quantities which are not involved in the back- 
ground configurations, as radial gauge field components 
and nondiagonal metric perturbations. It is found that 
this subset of perturbations is stable if the requirement 
of being bounded far from the hyperstring is relaxed (a 
condition which may not be required since there is not 
reference background fields for these perturbations). 

The article is organized as follows: after setting the 
field theoretic framework both for the particles and grav- 
ity in Sec. ^ we construct the Nielsen-Olesen ansatz 
for a three-dimensional vortex configuration, set and dis- 
cuss the corresponding Euler-Lagrange field equations in 
Sec, mil We then show how to handle the boundary con- 
ditions in Sec. IIVI and solve numerically the field equa- 
tions in Sec. insisting in particular on the numerous 
technical difficulties. This permits us to obtain the pa- 
rameter range over which gravity is localized and exempt 
of singularity in the core. We then discuss some argu- 
ments leading to the suggestion that such a defect real- 
ization of a brane in six dimensions may be marginal in 
Sec. lVIl and explicitly derive and discuss its allowed per- 
turbations in Sec. IVIll before ending by some concluding 
remarks. 



II. VORTEX CONFIGURATION IN ADSg 

We consider the action for a complex scalar field <& 
coupled to gravity in a six-dimensional space-time 



S ■ 



1 
2k 



2 (R - 2A) + Cm^t 



(1) 



where is the six-dimensional metric with signa- 

ture (-1-,—,—,—,—,—), R its Ricci scalar, A the six- 
dimensional cosmological constant and = 327r^Gg/3, 
Gg being the six-dimensional gravity constant^. The 
matter Lagrangian reads 

-Cmat = ^g^^T^A^ (2?B<i>)* - V{<i>) - \HasH^^, (2) 

in which capital Latin indices B . . . run from to 5, 
Hab is the electromagneticlike tensor defined by 



Hab — ^aCb QbCa 



(3) 



where Cb is the connection 1-form. The U(l) covariant 
derivative Va is defined by 



Qa - iqCA 



(4) 



where q is the charge. The potential of the scalar field 
<& is chosen to break the underlying U(l) symmetry and 



^ In D dimensions, we relate to the D-dimensional gravita- 
tional constant Go by k?^ = {D - 2)n'"-2lGD, where n'^-^' = 
IV{DI2) is the surface of the (D - l)-sphere. 
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thereby allow for topological vortex (cosmic stringlike) 
configurations, 



where a prime denotes differentiation with respect to r. 
Similarly, the matter stress-energy tensor. 



A 



(5) 



-L AB = r 4„ gABi-n 



(9) 



where A is a coupling constant and rj = (|<1>|) is the 
magnitude of the scalar field vacuum expectation values 
(VEV)^ 

Motivated by the brane picture, we choose the metric 
of the bulk space-time to be of the warped static form 

ds2 = g^^dx^dx^ = e''^'''>7]f,^dxt'dx'' - dr^ - r^e''^'''>d9^ , 

(6) 

where ry^j^ is the four-dimensional Minkowski metric of 
signature (-1-,—,—,—), and {r,9) the polar coordinates 
in the extra dimensions. Greek indices ^^v . . . run from 
to 3 and describe the brane world sheet and we set 
g^i, = exp[cr(r)]77p,y. The action ^ with the ansatz © 
will admit static solutions depending only on r so that the 
general covariance along the four-dimensional (physical) 
space-time is left unbroken. 

The Nielsen- Olesen like 64;| ansatz for a generalized 
vortex configuration is taken to be of the form [g^ 

$ = ^(r)e™^ = 77/(r)e™^ Co = -[n-Q{r)] (7) 

q 

where n is an integer, so that the only nonvanishing com- 
ponent of the electromagnetic tensor is Hgr = Q' /q- 
With such an ansatz, we shall now derive the relevant 
field equations and discuss their solutions. 



III. EQUATIONS OF MOTION 

With the metric given by Eq. lO, the nonvanishing 
Einstein tensor components reduce to 



has nonvanishing components provided by the Nielsen- 
Olesen ansatz {Tjl that are given by 



V - 
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2 j-2 



1 Q'2 



2 g2,.2 



2 T 

r e' 



-t; — e 



1 Q'' 



(10) 



It follows that the six-dimensional Einstein equations, 
with our conventions. 



Gab + J^9a 
can be cast in the form 

3 . . 1 .. 1 



^^JTab — Oj 



(11) 



3 .. 3.2 3 



1 



-7 = 



p W\ 

/?(/'-i)' + /' + ^[qV' + ^] 



,(12) 



3.2 2 . . . 

-a H — a + (7"f 
2 P 



A 



a 



(13) 



1„ ( I I c_'2 I Q„/„/ 



+27" -f 7'2 + -7' 
r 



\ ( 4 

G„ = cr' 3cr' + - + 27' 

2 V r 



Gfl 



1 



(8) 



^ Note, that, because of the unusual number of space-time dimen- 
sions, the fields have dimensions given by [ij] = M^, [<&] = M^, 
[A] = M2, [A] = A/-2, [r?] = Af2, [k^] = M"*, {Ca\ = M^, 
[Hab] = and [q\ = M^^ (M being a unit of mass). 

This can be further generalized in the D— dimensional case by: 
[$] = M(0-2)/2, [A] = M4--D, [,,] = M(D-2)/2^ [^2^] ^ M^-O, 

[Ca] = M°-2, [Hab] = M°-^ and [q] = M^-^/a^ [R] and [A] 
being unchanged. 



A 



/3(/'-l)' 



/ " — Q / + — 



(14) 



where we have introduced the dimcnsionless radial coor- 
dinate 



p=V\Mr, (15) 
as well as the dimensionless parameters 



1 2 2 



(3: 



1 A772 

4W 



? 2 



(16) 



In Eqs. (|12|l to p4|) . we have introduced the convention 
that a dot refers to a differentiation with respect to the 
dimensionless radial coordinate p. 



4 



The scalar field dynamics is given by the Klein-Gordon 
equation 

V^V-^* = -77^) $ + ^2^2^ (17) 

which takes the reduced form 
/ + (2a + + / = ^/c-'' + 2P if - 1) /, (18) 
while the Maxwell equations 

V^iJ-"^ = -g2c2|$|2 ^ * g ($a^$* _ $*9^$) , (19) 

provide the single reduced equation for the only nonvan- 
ishing component of the gauge vector field 

Q+(2&-^i-^^Q = efQ. (20) 

The set of equations (p H [T^ [HI [T^ is a set of 
five differential equations for 4 unknown functions {a, 7, 
/, Q). Indeed, we have a redundant equation due the 
Bianchi identities and one can check that the Higgs field 
equation (|18(l is recovered from the constraint equation 
(|13|l provided / 7^ 0. This system can be further sim- 
plified by remarking that Eq. H2()|l can also be written 
as 



d / e 



dp \^ym J i/to 



(21) 



Using the relation H24() with a = —b leads to 

i[e-VS(.-0]=i^(f + (28) 



and combining this result with Eq. (|23|l . one obtains 

4a gg 



s — I + c, 



(29) 



where c is a remaining integration constant |6a |. The 
equations of motion ijEl [El [H d HOJ) end up be- 
ing equivalent to the following five first-order differential 
equations 

5 ^ 2aQ'w ^1,; c 



4a ^ 

m = sm Qw -\- cm, 

£ 

So 4a Qw 



(31) 



V, V, 



af = ^s^-— — s + ^+^-^+c,5, (32) 
I e m 12 

w, (33) 

(34) 



2a Qw'^ 3 



w = sQf -sw + -w. 

e m 2 2 



where 



After discussing the behavior of these fields far from 
the vortex, i.e. far in the bulk, and on the brane itself 
in the following section, we shall solve numerically the 
field equations in order to determine the structure of the 
space-time and defect system. 



IV. ASYMPTOTIC BEHAVIORS 



so that 



2 2 'V 

m = p V = p e ' , 



( ^gg) =_(e/2g2 + Q2 
dp Vv"^ / v"^ 



(22) 



(23) 



On the other hand, combining Eqs. (|12|) and (|14ll . one 

also finds that 



dp 



(24) 



for any set of arbitrary constants a and 6, and where we 
have set the new functions 



as well as 



TO 2 

S = cr, t = — = h7, 

TO p 



T ^ a;3(/2-l)%A, 



V, 



2a g^ , PQ^ 
, V, = 2a . 

em TO 



(25) 

(26) 
(27) 



By definition of the topological defect like configura- 
tion, we require that the Higgs field vanishes on the mem- 
brane itself, i.e. <[> = for p = 0, while it recovers its 
VEV, 77, in the bulk. These requirements translate into 
the following boundary conditions for the function /: 



/(O) = 0, 



lim / = 1. 



(35) 



The corresponding boundary conditions for the 1-form 
connection are given by 

g(o) = n, g(o) = w(o) = 0, iim g = 0. 

(36) 

In order to avoid any curvature singularity on the 
string, the Ricci scalar stemming from Eq. (Q, namely 



i? = -|A| (7 + 4^ + ^72 



5a^ 



27 4(7 

270- + — + — 
P P 



(37) 

has to be finite at p = as the vortex is assumed to rep- 
resent our physical four-dimensional space. As a result, 
the warp functions 7(0) and (t(0) have to vanish in the 
string core. 



7(0) = d(0) = 0, 



(38) 
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and the warp function I therefore scales near the string 
hke 



l{p) ^ - 

p 



m oc p 




(39) 



Note that we impose both functions to vanish on the 
string, and not merely the combination j + 2a entering 
Eq. (|37l) : this arises from the requirement that all geo- 
metrical quantities, e.g. R^W^^ in which 7 and a enter 
with different coefficients, must be finite. A coordinate 
transformation along the brane allows to choose a{0) = 0, 
while 7(0) and Vg, defined by 



=7(0) 



(40) 



are determined by the boundary conditions at infinity. 
Note that cannot be absorbed in a rescaling of the 
radial coordinate. To see this, it suffices to introduce a 
new coordinate, f say, such that r^e'^^^^ = f^, which is 
equivalent to defining 7(f) — "f{r) — 7(0). This would 
induce a shift in the other warp function a, shift that 
can however be taken care of by a rescaling of the vor- 
tex internal coordinates. This is not all though, because 
the last metric element grr = —1 gets modified into 
gff = grrG'^'"^^ — — e''''-^-'. All the derivatives with re- 
spect to this new radial variable also acquire this numer- 
ical factor. In the Einstein tensor, given the symmetries 
of the vortex, this seems harmless as Gab is also sim- 
ply rescaled. However, the stress-energy tensor (|10|) is 
not so simply rescaled as it also involves non derivative 
terms (the Higgs field potential V and the gauge-Higgs 
coupling). On the other hand, Vg can be absorbed by a 
rescaling of the angular coordinate 9 — y/v^9. In that 
case, the angular part of the metric © appears to be 
cylindrical in the hyperstring core, with however a miss- 
ing angle 



A9 = 2^ (1 - ^) 



(41) 



The space-time geometry obtained for Vg ^ 1 exhibits 
a conical singularity in the vortex core whose physi- 
cal interpretation is the existence of an additional (5-like 
energy-momentum distribution (a Goto-Nambu hyper- 
string) lying at the center of the configuration. This in- 
terpretation remains valid provided < < 1, the other 
cases will be discussed in the next section. At this point, 
it is interesting to note that contrary to what is assumed 
in Ref. , the value of Vg in our approach is completely 
determined as soon as the other boundary conditions are 
imposed and ends up being a function of the model pa- 
rameters only. Setting Vg — 1 afterward, to obtain a 
regular geometry in the hyperstring core, will allow us to 
recover the fine-tuning relation obtained in Ref. |62| . 

In the following, we derive analytical approximations 
of the fields at infinity and in the hyperstring core asso- 
ciated with an anti-de Sitter space-time at infinity. The 
influence of the model parameters on these solutions is 
discussed. 



A. Far from the string 

Asymptotically, the anti-de Sitter space-time is recov- 
ered provided 



lim s 



lim I 



0. 



(42) 



Denoting by an index 'F' (standing for "fixed" ) the value 
of the fields at infinity, it follows from Eqs. H35|l . (|36(l and 
(|42|l that the adSg solution is a fixed point for the set of 
Eqs. 1201-123 with = 1, = = for the Higgs 
and gauge fields, and with the equations 



2 ^ 



£ _ A - n 
2*" |A| 



0, 



A 



2c.. + - = 0, 



(43) 
(44) 
(45) 



for the warp factors. The equations (|43|) and H45|) can 
only be simultaneously satisfied for c = 0, since Sp = 
would lead to A = 0. As a result, Eq. (|44() requires that 
TOp = and the asymptotic warp factors reduce to 



^" ^" 5|A|- 



(46) 



The anti-de Sitter solution is obtained for A < so that 
Zp = Sp = — 1/2/5 and we can now verify that and 
effectively vanish. Indeed, from Eq. H31|l . the dominant 
behavior of m at infinity is given by 



/2/5p 



(47) 



while Eqs. (|33|l and ()34|l admit asymptotically the decay- 
ing solutions 



with 




(48) 



(49) 



From Eq. (|?7|) . the gauge functions and vanish at 
infinity provided 



(50) 



This is the first restriction on the available parameter 
space. 

From Eq. H18|) . the decaying branch of the Higgs field 
at infinity is given by 



Arrir, 
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where we have defined 



h=l-f, 



and £h reads 




(52) 



(53) 



From Eq. (|51(l . it appears that the behavior of the Higgs 
field at infinity can be driven by the gauge field, provided 
£h > 2£g — a/2/5. In this case, one can check that h 
remains positive definite ensuring that the Higgs field 
approaches its VEV from below and thus can support 
a topological defect configuration. Indeed, would h be 
negative asymptotically, the condition /(O) = could 
only be realized if the sign of the derivative / changes at 
some intermediate point. From Eq. p8|l . it is clear that 
at the point / = 0^ the right-hand side is positive for 
/ > 1. As a result / can only be positive and the profile 
of the Higgs field would remain always convex and greater 
than its VEV. 

Similarly, the asymptotic expression for the warp fac- 
tor can be obtained from Eq. (|30|) . In terms of the new 
function 



u — s + 



(54) 



one gets for the decaying branch, up to some fine-tuning 
(see Sect. EH) 



Aaf3 



-2t 



24 



hP -( 2^,-^/275) p 
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2e TOc 



(55) 



If 2i?h < 2^g — •\/2/5, the convergence of the warp fac- 
tors toward the anti-de Sitter solution is driven by the 
Higgs field and u remains positive definite since hoc is 
positive in that case [see Eq. (|51|l ]. On the other hand, if 
2£h > 2^g — a/2/5, the metric factor s behaves asymptoti- 
cally like the gauge field and u remains definite negative. 
As the result, the surface 2£h = 2i'g — ■\/2/5, i.e. from 
Eqs. 63) and 



1 

10 



1 



7^' 



(56) 



separates the parameter space {a,e, (3) in two regions 
where the warp factor s approaches its anti-de Sitter 
value from above or below, respectively (see Fig.nj. 



B. Near the string 

The field behaviors near the string, to leading order in 
p, can be extracted from the equations of motion 13U|I to 
(|34|l . We assume the truncated Taylor expansions 



w 



(57) 



for the gauge fields, and 



(58) 



for the Higgs field with > 1 and k^, > 1 in order that 
w and / remains finite in p = 0. Similarly, the warp 
factor s can be expanded around p = as 



(59) 



with > 1 in order for the Ricci scalar to remains finite 
in the core [see Eq. ||2ZI)]. Setting c = (see Sect. II V A|l 
in Eq. using Eqs. gOl, (|SZI) and (jSHl, yields 

K — 1 and 



2a 



-nw. 



(60) 



From Eq. (|30|l and making use of Eqs. (|26|) and (|27|l . one 

gets 

At, = 1, s, = -{l-aP+ '-). 61 

The Higgs field behavior is given by Eq. 1321) and reads 

a«;^/>'"/-' = 2,s, + (a/3-1) 



? f9 ? 



With given by the expression H61|) . this implies 

n 



(62) 



(63) 



As a result, the behavior of the Higgs field around the 
string core is only determined by the asymptotic solutions 
at infinity through (and w^), as announced above, and 
the requirement > 1 only provides the constraint 



v„ < n 



(64) 



From a purely classical point of view, only / has to be 
well defined at p = 0. However, in the general situation, 
there always exists an integer p g N such that > /p, 
in which case all derivatives of the Higgs field, f^''\ with 
k > p are divergent in the core since the p^^ derivative al- 
ready is. This is so unless the bound is saturated, i.e. the 
equality — /p is strictly satisfied so that = p and 
then all derivatives of order larger than p will strictly 
vanish. If the bound is not saturated however, the met- 
ric at the core exhibits the conical singularity on which 
the field reacts by introducing the aforementioned diver- 
gences. Note that, as discussed above, the conical singu- 
larity interpretation demands < I, while for a vortex 
with 71 > 1, Eq. (|64|) allows to be larger than one. This 
unexpected regularity comes from the scaling properties 
of the field equations. Indeed, in Eqs. l|5n|) through |(33J, 
the warp function m (and thus v) appears only through 
the ratio Q^/m and w'^/m. As a result, a solution for a 
given winding number n and Vg , is also solution of a wind- 
ing number h = pn and = p'^Vg. This scaling permits 
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to solve the equations for a reduced set of parameters 
and yet obtain the complete spectrum of solutions. On 
the physical side, the kind of singularity appearing in the 
core for 1 < < may be interpreted as an supercrit- 
ical Goto-Nambu hyperstring lying in p = 0. On the 
other hand, the behavior of the fields leading to Vg > 
requires divergences of / and the stress tensor no longer 
remains finite in the vortex core [see Eq. (|10|l ]. In that 
case, the conditions > 1, > 1 and Eqs. are no 
longer valid and the regularity requirements in Eq. H38|) 
are no longer satisfied such that a curvature singularity 
appears in p = 0. 

The previous analysis makes clear that the value of 
Vg encodes the regularity of the matter fields and the 
geometry in the string core. By requiring the space-time 
to be of anti-de Sitter kind at infinity, we have shown 
that Vg is a function of the model parameters only, as 
long as the boundary conditions (|35|l . (|36|l and H4t)|l can 
be satisfied. In the following, we numerically recover the 
field behavior expected from the asymptotic analysis, in 
particular the regular configurations in the core require 
a fine-tuning in the model parameters such as = 1. 

V. NUMERICAL APPROACH AND PROBLEMS 

Several technical difhculties appear in the numerical 
integration of the equations of motion (|12|l to (|20|l . The 
first issue comes from the requirement of an anti-de Sit- 
ter space-time far from the core. Indeed, as mentioned 
in Sect. IIV Al the metric factors as well as the Higgs 
and gauge fields admit growing modes at infinity which 
correspond to infinite- volume space-time in the extra di- 
mensions (see Sect. I VII and Appendix IB)) . As the result, 
any direct numerical integration starting from the core 
toward the outer regions will necessary jump, due to the 
finite numerical accuracy, onto these growing solutions. 
This numerical instability can be overcome by perform- 
ing a backward integration starting from a finite cutoff 
distance far the hyperstring toward the core. However, 
one has to pay attention to choose a convenient set of 
functions, as s and m rather than a and 7 to remove 
any explicit 1/p dependencies in the equations of mo- 
tion. Indeed, the effect of the 1/p term in Eqs. (|12|) to 
(|20|l is to add flow turning points where the signs of the 
field derivatives change. As a result, growing behaviors 
would appear from these turning points whatever the ini- 
tial conditions and the direction of integration. There is 
no such flow inversions in the closed system of Eqs. H30() 
to H34|) . and the exponential growth can be suppressed 
by integrating from the anti-de Sitter fixed point at in- 
finity toward the hyperstring core. However, one has still 
to face the following difficulty. As pointed in Eq. (|55|1 . 
even the decaying branch of s at infinity admits two ex- 
ponential decaying modes: one varying as exp(— 2^hp) 
and the other as exp[—{2£g — y/2/5)p]. As a result, a 
backward direct integration is numerically unstable with 
respect to one of these modes. In the same way that a 



forward direct integration would jump onto the growing 
mode, such a backward numerical method would tend 
to select by numerical finite accuracy the strong est d e- 
caying exponential. For instance, if 2€h > 2^g — a/2/5, a 
backward method would be preferentially sensitive to the 
exp (—2ii^p) behavior (by moving toward the lower values 
of p, this mode blows up faster than the other), which 
is not the one in which we are interested. It is interest- 
ing to note that the numerical instabilit y of the back- 
ward method disappears for 2£h = 2£g — a/2/5 where the 
two decaying modes are identified, which means that this 
method can only be efficient on a surface of the parameter 
space {a,e,(3) whose equation is given by Eq. Note 
that a similar numerical instability occurs in Eq. H18|l be- 
tween the two corresponding decaying modes of Eq. H51|l . 

In order to overcome these difficulties, we have chosen 
to use a finite difference numerical method instead of a 
direct method. In particular, the equations of motion 
()f 2|l to (|2()|l deriving from the action we have im- 
plemented a successive over-relaxation method |67j . By 
discretizing the radial coordinate p, the action |^ can 
be expressed as a finite sum over pi (the integer i index- 
ing the discrete values taken by p over the radial grid) 
whose differentiation with respect to the fields evaluated 
at the discrete points leads to a system of finite difference 
equations corresponding to Eqs. ifT^ . (fTljl . ifTI^ and if^ . 
From a initial guess of the field and metric profiles along 
the radial grid, the values taken by the fields at pi are 
corrected by a Newton's method to reduce the error with 
respect to the true solution. In this approach, the bound- 
ary conditions are part of the finite difference equations 
since they appear as the conditions satisfied by the fields 
at the first and last point of the pi grid (see Appendix 
A in Ref. for a relaxation method applied to a simi- 
lar action). This procedure is stable provided the initial 
guessed profiles are not too far from the true solutions. 
The iterative corrections can been stopped when the dis- 
crete action remains stationary at the machine precision. 
We have also checked that the numerical solutions ob- 
tained in this way satisfy the constraint equation (|13l) . 

In order to probe the behavior of the solutions accord- 
ing to the model parameters, we have, in a first time, nu- 
merically allowed regular and conical solutions in the hy- 
perstring core by requiring the boundary conditions l|35(l . 
(ISnil, (ESi, llSni) and igni) to be satisfied. Along the lines 
drawn in the previous paragraph, we have used the relax- 
ation method to compute the solutions of the field equa- 
tions in the case of a unit winding vortex n — 1 and for 
25'^ values of the parameters {a,e,(3). As expected from 
the asymptotic analysis (see Sect. IIV|) the hyperstring 
generically develops a conical singularity with / 1 
(see Fig. ^ . The regular solutions Vg = 1 are obtained 
only for the parameters lying on the surface plotted in 
Fig. 131 which identifies to the fine-tuning surface previ- 
ously found in Ref. (62i] . 

However, under the previous boundary conditions, the 
relaxation procedure failed to converge in the regions of 
the parameter space which could have been associated 



FIG. 1: Typical field solutions associated with an anti-de Sitter space-time at infinity in the regimes where Vg < 1/4, 1/4 < 
Vo < 1 and Wo ~ 1, obtained for (a, e, (3) equals to (1.00, 5.00, 1.50), (2.20, 5.00, 0.80) and (1.29, 14.83, 6.39) respectively. In each 
case, the Higgs and gauge fields are plotted on the left picture, the warp factors are plotted in the middle one, while the right 
plot represents the behavior of f/p and v = m/p^. Note the behavior of the warp factor s (middle plots) for the parameters 
living above or under the plane 15fcill (see also Fig. Inland Fig.|HJ. 



with values of u^, > 1 (see Fig. For these regions, 
we have thus weakened the boundary conditions in order 
that Eqs. and (gBJ are stiU satisfied but with 



ct(0) = 0, m(0) = 0, (65) 



instead of Eqs. H38I) and H39|) . In that case, the method 
converges and the numerical solutions exhibit a diver- 
gence in the Higgs field derivative while the derivative of 
the warp factors do no longer vanish in p = (see Fig.0J. 
From Eq. I|37|l , these solutions correspond to a curvature 
singularity in the hyperstring core and are not physical. 



VI. FINE-TUNING AND STABILITY 



The solutions we have obtained for the fields surround- 
ing a branelike vortex requires a fine-tuning of the un- 
derlying microscopic parameters to be free of singularity 
in the string core and of finite volume in the extra di- 
mensions. These conditions are indeed the minimal re- 
quirements for an acceptable warped braneworld model 
in six-dimensions. In the following we discuss qualita- 
tively why such fine-tuning is expected as well as the dif- 
ferences appearing with respect to the domain wall model 
in a five-dimensional anti-de Sitter case p^ . Our analy- 
sis is analogous to the one developed in Ref. [s^, !54] in 
the case of global vortex and leads to similar conclusions. 
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FIG. 2: Isosurfaces of constant v„ in the parameter space {a,e,l3) associated with an anti-de Sitter space-time at infinity. From 
the left to the right, the four surfaces correspond to Vg — 0.25, Vg = 0.5, Vg = 0.75 and Vg = 1, respectively. The wired mesh 
represents the plane /3 = — l/10-|-e/4 (see Sect. 11 which separates the parameter space in two regions. For /3 > — l/lO-l-e/4 
the gauge field drives the warp factors toward their anti-de Sitter value and s has a global minimum at a finite distance to the 
core whereas, in the other case, the Ifiggs field dominates at infinity and s always decreases toward its asymptotic value (see 
Fig. 0. On the left side of the Vg — 1 isosurface the hyperstring exhibits a conical singularity in p = whereas on the right 
side there is a curvature singularity (see Fig. 2J . Only the surface Vg — 1 ends up be associated with a regular configuration on 
the brane (see Fig. |HJ 



The asymptotic form of Eq. H3U|) reads 



s + -s 

4 



1 

2' 



(66) 



for which one obtains the general solution in the form 



2 Ae^eP 

X -7^^ 

5 A.e^oP 



(67) 



where Cg = y^5/8, and is an arbitrary constant, to be 
matched with the vortex interior solution. In the asymp- 
totic analysis of Sect. IIVI this constant has been set to 
zero, but clearly, ii Ag ^ 0, one has 



lim s(p) 



(68) 



leading to an exponentially divergent warp factor for the 
metric ©. 

Only the particular value Ag = can smoothly join 
the interior metric to a six-dimensional anti-de Sitter 
asymptotic space-time. In five dimensions, this choice 
happens to be imposed by the Einstein constraint |23j |. 
but in the case at hand, because of the extra degree of 
freedom provided by the other function 7, this is an ex- 
plicit choice and not a mandatory consequence of the 
Einstein equations. In other words, the solution satisfy- 



ing hmp 



/2/5 is also a point in the phase space 



from which all trajectories diverge. This is related to the 
fact that in the limit where the field contribution in the 
stress-energy tensor is negligible with respect to the bulk 
cosmological constant, Eqs. (|12|l and H14|) are two dy- 
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namical equations for the warp functions cr and 7, with 
Eq. (|13|l being a constraint. The two first order (in a and 
7) dynamical equations require two constants of integra- 
tion, of which the constraint fixes only one (together with 
the requirement of an anti-de Sitter asymptotic space- 
time). The solution with arbitrary (nonvanishing) 
is thus a valid solution, contrary to the five-dimensional 
case. However, any nonzero value of does not corre- 
spond to a solution with gravity localized on the vortex, 
being exponentially far from the adSg case (see also the 
Appendix . One is thus led to conclude that the fine- 
tuning required in the 6D case is worse than in 5D since 
the physically relevant solution is a set of measure zero 
in the full set of solutions. This drives us to ask whether 
such a solution, although mathematically acceptable, can 
be reached by any dynamical evolution. In fact, as we 
show in the following section, no acceptable perturba- 
tion mode can be found for the regular vortex, except for 
some special modes, the transverse ones, which have no 
equivalent at zeroth order. In other words, the regular 
6D vortex configurations cannot be subject to perturba- 



tions in their background fields and in particular, cannot 
depart infinitesimally from adSe. 



VII. GAUGE-INVARIANT PERTURBATIONS 



As vector and tensor perturbations have been inves- 
tigated elsewhere 'gfl, |61|, we shall concentrate on the 
scalar part 63] of the perturbations. Here and in what 
follows, the Scalar- Vector-Tensor decomposition is un- 
derstood to be with respect to the four-dimensional vor- 
tex internal coordinates. Note that most of previous 
works were concentrating on zero modes, which were 
found to be "not normalizable" . In what follows, we 
consider massive modes and attention will be paid to the 
physical interpretation of the results. 



11 



0.5 



-0.5 











- 










s 




1,1, 


1,1, 



FIG. 4: Typical behavior of the Higgs field and metric factors 
for the model parameters living on the right-hand side of the 
= 1 isosurface, i.e. for large value of a (see Fig. 1^. The 
derivatives of the Higgs field and the metric coefficient diverge 
in p = leading to a curvature singularity [see Eq. 13711 ] . 



A. gauge-invariant variables 

In order to conclude on the stability (or physical rele- 
vance as we shall see) of the configuration we obtained, 
it is necessary to perturb this background solution in a 
gauge-invariant way. The first order perturbation of the 
metric ©, when restricting attention to the scalar per- 
turbations, reads 

ds2 = e-i-^) [^^, {l + i>) + df^d^E] dx^'dx" - (1 + e)dr2 

- 2Cdrde' - r^e^^'') (1 -|- uj) dO"^ 

- 2 {d^,Bdr + d^Cdd) dx", 

(69) 

where the scalar functions ip, E, ^, (, cu, B and C depend 
on all the coordinates (x" , r, 9) and are assumed to be 
small. 

A gauge transformation x^ x'^ — x'^ + e"*, with 
_ Qti^ (scalar transformations only) implies three 
gauge degrees of freedom, so we are left with four un- 
known functions to determine. The scalar functions 
transform under a gauge transformation as 



1 
2 



LU = UJ — 



- + 7 Cr, 

r 



B = B-^ (e, + ('- a'e) , (7 = C - i (ee + dge) . 

(70) 

From these relations, we can derive the four gauge- 



invariant variables 
^ = tlj~^a' {43 + e"E') , S ee ^ - 9r (4B + E') , 



T 



(:-\dr (4C + e'^deE) ~ ^dg (45 + e^E') 



-de (4C + e-'deE) 
7' ) (45 + e^drE) . 



^l — UJ - 
~ 2 

(71) 

As in the usual cosmological case pal: these variables 
are identical to the original variables once the choice of 
longitudinal gauge {E = B = C = 0) \s made. The 
transformation leading to this gauge, starting from an 
arbitrary gauge transformation, reads 

e = - ^e^E, e, = 2B + ^c^E', eg = 2C + ^e^dgE, 

(72) 

and is unique. This gauge choice, which we shall for 
now on adopt, is thus complete for metric perturbations, 
but also for the matter ones. Indeed, in our model ||2J), 
the matter perturbations concern only the hyperstring 
forming scalar field $ and its associated gauge field Ca- 
Note that in our framework, the location of the brane is 
given by the zeroes of the Higgs field and thus directly 
taken into account in its perturbations. Since we are only 
interested in scalar perturbations, the perturbed fields 
can be expanded as 



(73) 



where we have extracted the background winding phase 
ginS |.]-^g scalar field perturbations. Note that, for con- 
sistency, all the perturbations have to be invariant by 
a complete rotation around the hyperstring, and thus 
can be decomposed in discrete angular momentum modes 
around the string. Under the gauge transformation 
x"^ ^ x"^ — x"^ + e"** these perturbations transform to 



X = X - er<^ - m^r-egif, 



1? = ^? - ^Cgee, 



1 , 1' 



dr^^r- -^Cge'g + 2 — Cg - + i- £9, 



^-dg- CgSr ^Cgdgeg, 

[ . . (74) 

and are therefore not invariant. Similarly to the metric 
tensor decomposition, we then define the gauge-invariant 
quantities through the relations 



1 



me 



X^X-^^'i^B + e'^E') 



ip{4.C + e''dgE), 

(75) 
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for the Higgs perturbations and 

e = t?-i^Ce(4C + e"ae£^), (76) 

er = A- Y^Ce [4C' + {<fdeE)'\ 

+ + (4C + e-a9ii;), (77) 

Qe = 'de- ^C'e {4B + E') 

- Y^Ce {^deC + e^d^E) , (78) 

for the gauge field perturbations. They also match with 
the original variables in the longitudinal gauge [E = B = 
C = 0). 

Note that since we are interested in perturbation the- 
ory, we have to keep in mind that all the perturbed phys- 
ical quantities involved at some initial time have to be 



close to the background solution. This implies in partic- 
ular that we must impose on the physically meaningful 
perturbations to be initially bounded: of all the possible 
solutions of the perturbation equations which we discuss 
below, we shall retain only those for which neither long- 
nor short-distance, divergence appear. This, as it turns 
out, is extremely restrictive. 



B. Perturbed Einstein equations 

The Einstein equations, perturbed at first order, stem 
from Eq. (|ll|l . The perturbed metric tensor Sg^B is 
explicitly written in Eq. (|69|l and allows, by means of 
Eq. H73|) , the determination of the scalar part of the per- 
turbed Einstein and stress-energy tensors. They are de- 
rived in Appendix 1X1 and Eq. Hll|l leads, in terms of the 
gauge-invariant variables, to the following equations of 
motion 



n' - 3—a'deT 



+ 



3a" + 3cr'2 ^ 



1 7'\' /I 7' 



(*-s) 



e-T Q' 



dV 



+ <e"77^,<j - (9^ - dgSr) + ^'S' + —'PQ'S + —v^QdgA + —S 



j,2 q 

le-^ fq^ 
q 



2 \q 



for the (/i, i^) part of Eq. Hll|l. The (/i,r) and (/i,0) components lead to the following equations, respectively, 



(79) 



3^-' + n' 



-deT - -a 



3/1,7' 



1,17' 



r2 + 2K? 



^^{Qg-dee) + ^'S'} =0, (80) 



and 



dg (3* + S) - T' - a 



/ , 1 7' 



T + 2^2 1^ (9, - e') + <fiQA - q^^Qe^ = 0. 



(81) 



The purely bulk components of the first order perturbation of Eq. (|ll|l read 



-e""n (35- + r!) - 2^d^^ + 2^a'd9T- 



3a' + 2{- + t 



e-T Q' 



r 2 



dV 



<< - ^ — (9^ - deer) + f'S' - S - -^(pQdeA - —S 

\ r q r r dip 



le-"^ 



„ le-^ /Q'2 



2 r2 V g2 



^'Q' n + —qip^QQg ^ - AS = 0, (82) 
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for the (r, r) part, 



dip 



r2 - qip^QQg \ - r^c'^AQ = 0, (83) 



for the {9,8) component, while the mixed one (r, 6*) leads to the equation 



-e^'^DT + 2dg^'-(j'dgE 



qip^QQr i - AT = 0, 



dgm - I (4a" + 5a'2) T 



o~ — - <y 



1 



(84) 



r 



with "□" standing for the flat four-dimensional 
d'Alembertian, i.e. 



(85) 



and where the perturbed Higgs field has been decom- 
posed as 



X = j:^iA, 



(86) 



i.e. into its real and imaginary parts. 

The perturbation JHEJ, although true in general, does 
not make the most of the U(l) invariance of the the- 
ory Indeed, with the definition Q for the covari- 
ant derivative, the full theory is unchanged under the 
change <i> ^'^^ = e'"'^'*'^ provided the gauge vec- 
tor field is simultaneously modified into Ca ci^-* — 
Ca — (l/q)dAa- For an infinitesimal U(l) gauge trans- 
formation, this leads to the transformation Z'^'^' ~ S 
and Z\(^^ = A + aip{r), so that choosing a = —A/ip al- 
lows to restrict attention to purely real perturbations of 



the scalar field perturbation X; we shall accordingly call 
this choice the real gauge. Indeed, as we shall see ex- 
plicitly, letting A arbitrary leads to equations involving 
only the U(l) gauge- invariant degrees of freedom, namely 
Q-A/{qip), Qe~deA/{q^) and e,-(l/g)(zi/^)', which 
merely expresses the fact that by going to the real gauge, 
one can get rid of A. 



C. Perturbed Maxwell equations 

By means of Eq. (|73|l . we can also derive the per- 
turbed Maxwell equations stemming from Eq. H19|l . at 
first order in the fields. Since the Einstein equations im- 
pose to the stress-energy tensor to be conserved, the per- 
turbed Maxwell equations are certainly already included 
in Eqs. I|79(l to (|84|l . Nevertheless, they mainly involve 
the matter fields and may help to decouple the whole 
system. The (/i) component of the perturbed Faraday 
tensor gives 



e; - e" + ( a' + i + ^) (e. - e') + ^-^dg (e^ - a^e) - q^A -t- gV'© = o. 



(87) 



while the (r) and (ff) bulk parts lead to 



e-'^n(e,-e') + _ae(e^-a,eo 



e-^ Q' . 

~ ' 

q 



2*--(S + 17) 



-q^'^QT - qipA' + qip' A 
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and 



Q" 



(^ 


4)] 







(e^ - deQr) + — 
q 



— S - 29^Qi: - qipdeA + gV'0e - 0, 
9 



(89) 



where use has been made of Eq. H2U|) to simplify the term otherwise proportional to f2. 



D. Perturbed Klein-Gordon equation 



In the same way, the Klein-Gordon equation l|17|l can also be perturbed in terms of metric and matter fields. By 
means of Eqs. ()73(l and JHSJ, its real and imaginary parts lead to two coupled equations, 



ip' 2*' + 



1 r 

r 2 
17' -S' 



E" -[2<j' + - + - ^—dlE + 1—QdeA 



A 



E 



-deT 



2a' 



1 , 7' 



^^Q^n - 2^q^Qee = 0, (90) 



and. 



1 , r 



e-^DA - A" - {2a' + - + ^] A' - ^d^gA - 2^QdgE 



+ (pQ- 



r 2 
T' -del 2^) 



r 



A 



qip 



+ 



ipQ' + 2^'Q 



9a' ^ ^' 
Za — 



-T 



V+(2a' + i + ^)^ 



qQr = 0, 



(91) 



respectively. 



As previously noted, the perturbed fields and geometry 
have to be invariant by a complete rotation around the 
string, i.e. they are 27r-periodic in the angular variable 6. 
Therefore, they can be decomposed in Fourier series with 
respect to 9, e.g. the perturbed Higgs fields is expanded 
as 

Z\(r,0)=^Z\p(r)e^P^ 

(92) 

i;(r,0)-^i:p(r)e*f^ 

and similarly for all the other perturbations. Plugging 
Eq. (|92|l . and analogous mode expansion for "if, S, fl, 
Q, Qr and Qg into Eqs. (|7^ to clearly shows that 
each angular "p-mode" decouples from the others. As a 
result, the time evolution of the perturbations can be fo- 
cused on a particular angular mode p, the physical evolu- 
tion ending up be given by their superposition. Although 
some redefinitions of the fields may separate the system 
of equations (|75|) to l|nT|l into distinct subsets 0, such 
a situation already happens for the lowest angular mode 
p = 0. Indeed, the angular dependency of these zero 
modes disappear which is formally equivalent to nullify 
the operator in Eqs. GHJl to One gets two 

disjoint pieces of equations, namely Eqs. (|79|l . (|80|l . (|82|l . 



(|5S|l . (|5^ and lEOI) which only involve Q, S, Qg and 
E on one side, while Eqs. l|HT|) . l|g7|l . (jS5|) and ^ 

couple only T, 0,., Q and A on the other side. 

In fact the zero angular momentum modes, obtained 
for p = 0, represent cylindrical perturbations which 
strictly wind around the string as the background form- 
ing fields do. Moreover, since they correspond to the 
lowest angular momentum state, one may naturally ex- 
pect them to be first excited in a generic modification of 
the vortex structure. 

In order to study the stability of these zero-modes, the 
time evolution of the latter subset of perturbations will 
be thoroughly analyzed in the following section. 



E. Stability of the transverse perturbations 

The time evolution of the lowest angular momentum 
modes , 0^ , and A^ is readily governed by the per- 
turbed Einstein equations H81I) and (|84|) together with the 
perturbed Maxwell equations H87I) and H88|l . and the per- 
turbed Higgs one As noted before, due to imphcit 
stress-energy tensor conservation in the Einstein equa- 
tions, this system is not over-determined although it in- 
volves redundant equations. Moreover, since we are inter- 
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ested in perturbations which behave almost hke the back- 
ground vortex fields, we will only consider real pertur- 
bations of the hyperstring forming Higgs field, i.e. with 
— (in other words, we go to the real gauge). In 
terms of the dimensionless background fields and param- 
eters [see Eqs. (0, (fT5|l and (fT?^ ]. together with the new 
dimensionless fields 



^ ge, 0r = 



q 



Or 



T= V|A|T, 



(93) 



the time evolution equations stemming from the Einstein 
and Maxwell equations read 



r - 



T 



4aQ 



(Or - ^) + 4:afQe 
f - AafQOr 



Or-0 



0r-O 



M-" Or-0] -S 



T 



while the Higgs one becomes 



m 



Q + Q\2 



: + 2s- 




0,(94) 
0,(95) 
0,(96) 
0,(97) 



T 

m 

(98) 



A four-dimensional Fourier transform has been per- 
formed on the zero angular momentum perturbed fields 
with respect to the four-dimensional coordinates x'^ . The 
general perturbation solution is therefore a linear super- 
position of the d'Alembertian eigenmodes 



r„{x'',r)= / T(/c^,r)e-*'=^^ d^fc 



(99) 



with f{kf,r) the solution of Eqs. ^ to and 
similarly for the others perturbed quantities. The 
d'Alembertian eigenvalues ends up being 

□ — > ~\A\M^ = -fi^'^k^K, (100) 

and any perturbation with positive mass squared > 
will be considered stable, whereas tachyonic modes, 
having Af^ < 0, will generate instabilities. 

There are three variables Or, O and T for five equa- 
tions, two of them being thus constraint equations. By 
means of Eq. (|95|l . the metric perturbation T can be ex- 
pressed in terms of Or only 



T = 



AapQ 



Or, 



(101) 



while by means of Eq. H101|) , Eq. H97fl gives the relation 
e^ = (p + l)0, (102) 



with 



2n2 



(103) 



Finally, plugging the previous expressions for T and Or, 
given by Eq. (|101|l and Eq. H102(l . into Eq. (j^El), one gets 
a second order differential equation involving only the 
function O, namely 



VO 



V 



V 



+ efO^Q. 



(104) 



One can also verify that the two remaining equations 
and also lead to Eq. (|104|l when use is made 
of Eqs. H101() and l|102|l . ensuring the consistency of the 
gauge choice — 0. If this choice is relaxed, one can 
check that the new perturbed equations simply require 
O to be replaced hy O— A^/ip, and Or by Or — dplA^/ip) 
in Eqs. (|1U1|I to (|103|l . The additional perturbed equa- 
tion stemming from Eq. (|91|l ends up be equivalent to 
Eq. (|104|1 . As expected for a gauge degree of freedom, 
the field A^ has therefore no dynamics and will not be 
considered in the following. 

In order to conclude on the stability of the vortex solu- 
tion with respect to these perturbations, let us consider 
the generic case of a real squared mass G M. In this 
case, and far from the string, Eq. (|104|l behaves as 



Af^ exp 




(105) 



and through the change of variable z = exp(p/vlO) and 
function O = exp(-\/5/8/9)'7', Eq. H105(l reads 




lOAf^ 



25 

4^ 



r = o. 



whose solutions are known, see, e.g. Eq. 
Ref. [tiJ. This gives 



O (X exp 




X Zi 



5/2 



10|Ar|exp 



(106) 



.491) in 



(107) 



in which Z5/2 is a Bessel function of order 5/2 and of 

its argument in brackets for Af^ > 0, and a modified 
Bessel function for M"^ < 0. As a result, for any positive 
mass squared, Af^ > 0, the solution given in Eq. H107() 
behaves, asymptotically far from the vortex, as an oscil- 
latory exponentially divergent quantity whose amplitude 
scales as 



\0\ oc exp I \j -p 




(108) 



Such solutions appear to be unbounded far from the 
string and are, strictly speaking, not well-defined: the 
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first order perturbation equations can be derived from 
the action expanded to second order in these perturba- 
tions, which one would thus expect to be finite, since 
the volume of the extra dimensions is finite (this is the 
very reason for choosing anti-dc Sitter in the first place). 
But this action contains a term cx / —geT" (the ex- 
ponential stemming from the unperturbed contravariant 
metric coefficient) which diverges exponentially. In spite 
of this issue, the solutions (|107|l can be given a physical 
meaning in the framework of perturbation theory since 
the energy they contain to first order is actually finite, 
and the corresponding gravitational potential T vanishes 
asymptotically [see Eq. H101() ]. In this respect, these posi- 
tive squared mass perturbations are physically admissible 
solutions. 

In the case < 0, the Bessel function in Eq. H1U7|) 
is of the modified kind and admits an asymptotically ex- 
ponential of exponential decaying behavior: one of the 
two degrees of freedom of Eq. H1U7|) has to be fixed to en- 
sure the decrease of the < solution asymptotically. 
However, this does not mean that tachyonic modes exist 
inside the system since it is also necessary that these per- 
turbations are well-defined in the string core. We shall 
accordingly turn attention to the interior solution. 

In the hyperstring core, the function V can be ex- 
panded as 



(109) 



where use has been made of the background field behav- 
iors in the string core. In the limit p 0, Eq. H1U4|) 
becomes 



(110) 



whose solutions are 



-Zi ( ^|M2-4a/2|p 



(111) 



Again, Zi refers to the two independent Bessel functions 
Ji and Yi provided — 4a/^^ > 0, while it designs the 
two modified one, Ii and Ki, in the other case (which is 
the case for AP < in which we are interested). There- 
fore, there always is a divergent solution in the string 
core, behaving as 1/p^, together with a well-defined one 
going to a constant value. As a result, the NP < de- 
creasing solution far from the string can not generically 
match with the well-defined one in the string core, the 
required degree of freedom being already fixed to ensure 
the asymptotic normalizability. Although there is thus 
no tachyonic continuum spectrum, the matching between 
the two well-defined solutions at infinity and in the string 
core could happen for some peculiar values of the nega- 
tive mass squared, making a discrete spectrum of unsta- 
ble modes. In the following, we show that this is not the 
case. 



The equation (|104|l can be rewritten in the form of 
a zero mode Schrodinger equation for the quantity u = 
[plVl exp (cr -I- 7/2)]^^^ 0, namely 

-u + Vm{p)u^O, (112) 

with the potential 

Vm{p) =W^ + W + ef + 4.a^—^ - M^e-", (113) 



with 



W=- 



m 



(114) 



the superpotential-like of the Schrodinger equa- 
tion (Ell "zl. 

It is immediately clear from Eq. H113(l that the poten- 
tial is asymptotically dominated by the last, exponen- 
tially increasing, term. As a result, a confinement for 
the corresponding scalar mode is achieved provided this 
last term is positive, hence requiring a negative squared 
mass. This would seem to imply the existence of tachy- 
onic modes on the hyperstring. However, a closer ex- 
amination of the potential (displayed on Fig.jS)) actually 
shows that even in the negative squared mass case, Vm is, 
numerically, positive definite: the associated Schrodinger 
equation has only strictly positive eigenvalues, and in 
particular no zero mode. This can be seen analytically 
in the following way: we first note that Eq. H112|l can 
be written as A + Z^{p) u = 0, where Z can be ob- 
tained from the potential H113|) and the definition of the 
operator 



A 



dp 



dp 



W{p), (115) 



APe' 



which is indeed a positive definite function of the core 
distance. Then, if one demands the perturbation u to 
be bounded, as implied by the fact that the energy con- 
tained in the mode be finite, then one is seeking zero 
mode bound (normalizable) states of the Schrodinger 
equation (|112|l . It is however immediately clear that since 
the operator A'' A has only nonnegative eigenvalues, then 
the spectrum of A^A + Z^{p) must be positive definite; 
thus, as announced, there are no zero modes solution of 
Eq. H112() . and hence no instability. It is interesting to 
notice that the nonnormalizability property at the hyper- 
string core also stems from the remark that the potential 
(|113|l diverges as Vm ~ near the core. 

Thus, the tachyonic solutions cannot be considered as 
physical perturbations of the system. We are led to the 
conclusion that, for the subset of fields considered in this 
section, the hyperstring is stable with respect to trans- 
verse perturbations. To decide on the overall stability of 
the vortex, one needs to clarify the role of the comple- 
mentary subset of modes; This is done in the following 
section. 
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FIG. 5: Characteristic shape of the potential Vm{p) given 
by Eq. 111311 for a typical set of values {a,l3,e) — 
(1.77, 2.66, 11.56) on the fine-tuning surface for which i;,, = 1 
and different values of the negative squared mass . It is 
clear from this figure that the potential is positive definite, 
and its minimum increases with \M^\. Hence, there is no zero 
mode for the perturbation equation, and therefore no insta- 
bility zone. 



F. Sausage perturbation modes 

We now turn to the second decoupled subset of cylin- 
drical perturbation modes which respect the axial sym- 
metry, hence their name. Provided we are interested in 
the zero angular momentum modes, the time evolution 
of the perturbations S, Q, E and Qg is determined by 
the Eqs. ^ and Some rapid 

simplifications can be performed. First, from the ^ ^ v 
part of the Einstein equation H79() . one gets 



(116) 



This expression can thus be used to simplify the above 
mentioned equations for fi — v. In terms of dimensionless 
quantities the perturbed Einstein equations simplify into 
two dynamical equations 



(2V, + 2VJ * + (2.F + V, + 2V,) S = S^^, 
10s* - 2sE + (Se'-^M^ - 2V, - 2V,) ^' 



4* 



and two constraint equations 

*-s + (.5-0^'-(s + /)s = V- 

2(s -I- l)i> - 2sk -I- {&-'^M'^ - 2Vi + 2V, 
-c-'^M'^ + 2J^ - V, -t- 2V, 



(117) 



(118) 



(119) 



(120) 



where the matter source field functions S stand for 



= 4a^©e - 4a//^ + 2(V, + Vj^e 
em 



Aaf + 2V, + 2/ 



df 



(121) 



-4a— Oe - 4.affIJ ~ 2(V, + Vj0e 
em 



Aaf - 2V, + 2/ 



d/ 



-Ob ~ 4ajjS, 



em 



(122) 



(123) 



em 

A T 

-4a/2 + 2V, + 2/— S, 



(124) 



in which 0q = qQg/Q and S = I^/ip. 

The perturbed Klein-Gordon and Maxwell equations, 
also expressed in terms of dimensionless quantities read 




+ /3/(,r-i) 



+ - Og) , 

m \ / 



(125) 



and 



Og 



2^ + 25-lj Og + e~^M^0g 



3^*-e/2 (e + 2S 



(126) 



It is interesting to note at this point that the equations of 
motion for the five-dimensional domain wall having this 
kind of perturbation modes are directly obtained from 
Eqs. ifTTTjl to ifT^ by setting Q = and removing any 
dependencies in I and Q. 

As previously noted, due to the Bianchi identities some 
of these equations are redundant. Indeed, Eq. (|117|1 is 
readily obtained by differentiation of Eq. (|119|l . up to the 
background Einstein, Klein-Gordon and Maxwell equa- 
tions ijSni) to (PHl. Similarly, differentiating Eq. (|120|1 . 
and using Eqs. (|119|l . (|125|l and H126I) to express Srr 
in terms of Sgg and the metric perturbations leads to 
Eq. H118|) . also up to the background Einstein, Klein- 
Gordon and Maxwell equations, and provided 1^0. As 
a resuh, only Eqs. fTT^ . fTB|) . ^2^, and 
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are relevant for this subset of matter and metric pertur- obtained, as it was the case for the transverse modes [see 

bations, with the constraint that the solutions are regu- Eq. I|l()4|l ]. 

lar at the point where I vanishes. However this system To study the stability properties of the hyperstring 

remains fully coupled and no simple second order dif- with respect to the sausage modes we derive in the foUow- 

ferential equation on one perturbation variable can be ing their behavior in the string core and asymptotically. 

I 



The previous system can be recast into a set of first order differential equations, 



21E = 



-e-^M"^ + - 2V., + 2s(s - I) 



/ 



AT 



4aff \^ + 2s\-2V,- 2f— 



f 



d/ 



S + iaffX + [ - 2V, + 2V^ - 2s{s - l)]0 - {s ~ l)W, 



-e'^M"^ + 2V, - 2V, + 2(s + l){s - I) 



e'^M^ + V,- 2V., -2T - 2{s + if 



/ 



2s + 21] -2V^~ 2/ 



df 



(127) 



S + iaffX +[-2V,+ 2V, - 2{s + 0(s - l)]e - {s - l)W, (128) 



for the metric perturbations 

X = 



2— + 

m f 



2— - ^(s + l) + — — 



f 1 
2- + 2.S + -I 
f 2 



X 



m / 
m / 



2a d/ 



2^ + Us -I) 



0, 



e-'^M^ + Aaf + 4/3/' 



(129) 



(130) 



for the Higgs field perturbations and 



e = w, 
.Q 



21W = 3 



Q 



e-"M^ + 2Vi - 2V, + 2s{s - I) 



-2elf 



Q 



.4d/2 + 3§ 



/ 



iaff l^ + 2s\-2V,- 2f— 



f 



df 



S+12a^ffX 



216-" Ar + 3^ [-2V, + 2V^ - 2s{s -l)]>0 



Q 



(131) 



e-^K'r + V, - 2V, -2T~ 2s{s + I) 



{3s + I) + 2l{2s - -I) 



W, 



(132) 



for the gauge field. 



From Eqs. (|127|1 . (|128|1 and making use of the asymp- 
totic behaviors of the background fields (see Sect. lIVp . 
the metric sausage perturbation modes at infinity verify 



and 



-eV^^M^- 




2 2 



(134) 



From Eq. H133|l , one gets 



/5''M2 + 2/5 



while Eq. (| 134(1 yields 



(135) 



(136) 



From Eq. p()5|l and ((l()7|l . the metric perturbations 
and therefore S, diverge at infinity as exp {yj2/hp) as 
long as M"^ > 0. This behavior is not admissible in 



19 



the framework of perturbation theory. As can be seen 
in Eq. (|69|l . the metric perturbations \I' and S have to 
be smah compared with their corresponding background 
values, themselves are of order unity, at least initially. 
Otherwise, mathematically speaking, it is not consistent 
to expand the equations of motion to first order. Phys- 
ically, this means we would start from a space which is 
infinitely far from the background one. 

Similar conclusions also hold for the matter fields: from 
Eqs. (|129|l to (|132|l . one gets asymptotically for the Higgs 
field perturbations 

S - + eV^PM^S ~ 0, (137) 

and, using Eq. 




V5P]^2^ ~ -3^g*, (138) 



for the gauge field. For AP > 0, both of these equ ation s 
have only divergent behavior at infinity, as exp ( y/2/5p) 
for the Higgs perturbations and as exp [{2ig + y/2/5)p/2] 
for the gauge field perturbations [see Eq. H5()(l ]. Also for 
the matter fields, the > solutions are not admissible 
since it would physically means that the Higgs field is 
infinitely far from its vacuum expectation value / = 1. 

As a result, the asymptotic study of the sausage modes 
shows that the hyperstring cannot be stable with respect 
to these perturbations. As this stage, either there are 
instabilities if there exist some < modes which 
are well-defined in the hyperstring core and match the 
decreasing solution at infinity (see Sect. IVII E\i . or the 
configuration is not perturbable at all, i.e. the only ac- 
ceptable solution is vanishing perturbations. To explore 
this point, we discuss the behavior of the solutions in the 
hyperstring core. 



J 



The physical solutions we are interested in have to be 
well-defined at p = 0. In particular, the geometry can 
only be regular provided ^'(0) = S(0) = [see Eq. l|57|) 
and discussion below for the background case] . Moreover, 
the sausage perturbations are required to be small with 
respect to their corresponding background values, i.e. 
S, S and O have to be finite in the core. Assuming these 
fields can be expanded in Laurent series around p = 0, 
the previous constraints yield 



n=2 



n=2 



(139) 



n=0 



n=0 



where ■(/)„, (t„ and 6'„ are real numbers. Plugging these 
expansions into Eqs. ((T77jl . ifHUI) . lfnT|) and 

H132|l . where the derivatives with respect to p are readily 
obtained from Eq. 1)1391) . leads to a set of coupled alge- 
braic relations for the coefRcients V'm ^n, c„ and On- By 
using the behaviors of the background fields in the hy- 
perstring core obtained in Sect. II V HI we find that this 
hierarchy is fully determined by the knowledge of the 
three parameters V'Oj ctq and 02- As a result, the reg- 
ular solutions in the hyperstring core generate a three- 
dimensional subspace of the six-dimensional full space of 
solutions. It is therefore necessary to fix three degrees of 
freedom to get regular solutions in p = Q. 

From Eqs. l(T^ . ((nS|) . ((nTjl and lfnS|l . we see that 
three degrees of freedom must also be fixed to ensure 
that the sausage perturbations are asymptotically well- 
defined (one for the metric perturbations ^ and S, one 
for the Higgs field perturbation S and one for the gauge 
field perturbation O). Moreover, from Eqs. (|127|l . H128() 
and H132|) . there will be no jump in the derivative of the 
perturbations at / = provided 



(-e-^M^ -I- 2V, - 2V, -f 2s2) ^,{p^) + (^-e-^M^ + - 2V, -2T- 25^) S(pc) 

^(pc) + iaffXipc) + [-2V, + 2V, - 2s2] 0{p,) - sWipc) = 0, 
I 



4affU + 2s]-2V,-2f^ 



(140) 



where the background fields are evaluated at p — pc, 
vanishing point of I. 



the 



Three degrees of freedom are thus fixed to ensure a 
convergent behavior of the perturbations at infinity, plus 
another one for the regularity at ^ = 0. Two degrees 
of freedom are left, which is not sufhcient, according to 
the previous discussion, to ensure regularity of the solu- 
tions in the hyperstring core. Even if one tunes the mass 
to keep only the solutions regular at p = 0, the con- 



vergent solutions at infinity (and regular at Z = 0) will 
not generically match with the regular ones in the vortex 
core. We have also numerically verified that no excep- 
tional hidden symmetry realizes this matching for a wide 
range of masses. Nevertheless, note that another param- 
eters of the model could be used to realize the matching 
between the regular solutions in the core and at infinity. 
Indeed, the background fine-tuning between a, (3 and e 
(see Fig. O is a surface in the three-dimensional parame- 
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ters space and one cannot exclude that instabilities could 
marginally occur on a curve along this surface. 

In conclusion, the only acceptable sausage perturba- 
tions modes are the vanishing ones, i.e. the hyperstring 
cannot be perturbed at all in the subset of matter and 
metric perturbations which correspond to nonvanishing 
background fields. In the following, we discuss the phys- 
ical meaning of this result. 

From a geometric point of view, it is well known 
that the generic space-time generated by a cosmic string 
in presence of a cosmological constant is of infinite- 
volume [7^ (see Appendix^ for a six-dimensional anal- 
ogous derivation). This is precisely why we have to 
fine tune the model parameters to obtain a finite- volume 
space-time with decreasing warp factors and no singular- 
ity in the core. As can be seen from the metric (0, the 
obtained space-time geometry leads to vanishing proper 
length circles around the hyperstring at infinity. This 
kind of geometry implies the existence of a point where 
l{pc) = which is precisely the stationary point of proper 
length circles. From p < pc the proper perimeter of a cir- 
cle around the hyperstring increases with respect to the 
radius, whereas for p > pc it decreases toward (see 
Fig. To link this structure to the usual conical ge- 
ometry generated by cosmic strings, one may imagine a 
missing angle starting from zero in the hyperstring core 
toward 2tt at infinity. This is not really surprising since 
we have required the hyperstring to generate an anti-de 
Sitter space-time at infinity, or naively, the fine-tuning 
allows to pass from a cylindrical symmetry in the core 
to a spherical one asymptotically. Now, it is clear that 
disturbing the fields around the values which lead to such 
fine-tuned gravitational configuration is not necessary al- 
lowed. And this is precisely our result. The only allowed 
perturbations concern the transverse modes which have 
not equivalent at zero order. On the contrary, all pertur- 
bations of the background fields, those generating this 
fine-tuned space-time, are forbidden. Interestingly, the 
natural behavior of the perturbations far the string for 
> 0, as exp(-y/2/5p), looks like the generic metric 
coefficients which appears when there is no fine-tuning 
(see Appendix IB|| . Although one might expect the hy- 
perstring to relax toward this generic configuration, no 
conclusion can be drawn from the perturbations theory 
since the generic space-time configuration, with infinite- 
volume, is not at all "close to" the studied fine-tuned 
one. To end this section, it is worth pointing out that 
the previous conclusion is valid in the physical motivated 
framework where the vortex remains regular in p — 0. 
This hypothesis has allowed us to set the regular expan- 
sions in Eq. H139() . On the other hand, as shown for the 
background fields (see Sect IIVI and Sect.0, there is a 
dense set of solutions associated with an anti-de Sitter 
space-time at infinity which exhibit conical, or curvature 
singularities, in the vortex core (see Fig. 0J. This sug- 
gest that an allowed evolution of the fine-tuned regular 
vortex might be the birth of a singularity in p = (see 
Appendix ^ . 



VIII. CONCLUSIONS 

In the braneworld framework, one issue is to determine 
how to model the brane, and in particular to investigate 
whether its internal structure influences the properties of 
gravity and of the other fields living on the brane. Among 
other solutions, more interests has been focused on the 
possibilitvfor the brane to be realized by a topological 
defect 0,11 El • 

In five dimensions, it has been found that there al- 
ways exists a domain wall solution that confines grav- 
ity, which moreover is symmetric with respect to both 
sides of the brane provided the usual relationship be- 
tween the bulk and the brane cosmological constants is 
satisfied. This relationship translates into a fine-tuning 
of the underlying microphysics parameters [23!|. As far 
as the gravitational sector is concerned, the properties 
of the braneworld are mostly independent of the inter- 
nal structure of the brane. It is then possible to find 
various confinement mechanisms that lead to the exis- 
tence ofbosonic P^illlilllliiliilsilslllllll as 
well as fermionic zero modes that can be made mas- 
sive jH (although with a spectrum not yet compat- 
ible with accelerator data). In short, a five dimensional 
topological model of our Universe is, for the time being, 
an open possibility both from the cosmological and par- 
ticles physics points of view. 

In six dimensions, the general machinery used to study 
five-dimensional reflection symmetric braneworld does 
not apply, mainly because of the necessity to regular- 
ize the long range gravitational self- interaction [43. A 
way around is to specify a complete model determin- 
ing the internal structure of the brane in order to grasp 
some features of six-dimensional braneworld models. In- 
deed, one will then need to discuss the genericness of the 
conclusions drawn on a particular microphysics. For in- 
stance, there exists a vortexlike brane configuration on 
which gravity was shown to be localizable Jg^, . We 
have shown in this article that such vortexlike branes are 
generically associated with a singularity in the core, ex- 
cept when a fine-tuning between the model parameters 
is assumed 62] . 

In order to address the stability issue of the fine-tuned 
solution, we have performed a full gauge-invariant per- 
turbation theory around the regular six-dimensional vor- 
tex background solution. Focusing on the scalar pertur- 
bations, we showed that the hyperstring forming fields 
and nonvanishing metric parts cannot be perturbed at 
all. This result comes from the requirement that the 
hyperstring generates a finite volume space-time with in- 
finite extra dimension and remains regular in the core. 
Any perturbation of the background fields would de- 
stroy this configuration and are not allowed. As a re- 
sult, a nonempty universe, where additional observable 
fields would generate such perturbations, is severely con- 
strained, if not altogether ruled out. Indeed all induced 
modifications of the string and gauge forming fields, as 
time and radial metric factors are forbidden, even at the 
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perturbation level. In other words, such nonempty uni- 
verse should not couple to Higgs and gauge fields, and 
should not modify at all the gravity in the radial extra 
dimension. Therefore, the physical status of such a con- 
figuration seems rather unclear, rather artificial to say 
the least, and the possibility of having such a 6D vortex- 
brane realization in nature very dubious. 

A priori, these conclusions are specific to the case at 
hand. In particular they might be argued to depend 
on the field content of the underlying theory. However, 
since the discussion involved the gravitational sector, it 
could be conjectured that six-dimensional braneworld 
with anti-de Sitter bulk and infinite extra dimension, can- 
not similarly be perturbed without exhibiting singulari- 
ties. This would imply that nonempty multi-dimensional 
braneworld models could only have one (large) extra di- 
mension of the warped form (unless some extra structure, 
such as a 4-brane at a fixed and finite location away from 
the 3-brane is added Q)- 
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By means of 

^5"" = (A2) 
one can get the inverse perturbed metric tensor 
5g^"' = -e-''r]^"'^, Sg^''' = 0, Sg^''^ = 0, 

p-7 p-7 (A3) 

5g^'^'—T, Sg^^ = E, 5g'' ^"—n. 

The perturbed Ricmann tensor can also be expressed as 
a function of the perturbed metric tensor through the 
perturbed Christoffel symbols 

5R\^„ = -5V\^.^, + ^r-^,,^^, (A4) 
where the covariant derivatives with respect to the un- 
perturbed metric have been noted with a semicolon, and 
the perturbed connections are given by 

S^\c = ^5-^" {^9oB,c + 6goc-B - Sgsc;o) ■ (A5) 



2. Einstein tensor 



APPENDIX A: PERTURBED QUANTITIES 

In this appendix, we derive all the gauge-invariant 
parts of the various tensors necessary for the stability 
analysis. 

1. Metric tensor 

According to Eq. Ht)9|l , the scalar perturbed metric ten- 
sor reads, in terms of gauge-invariant variables, 

'^ffM'' e'^jy^^^', %r = 0, %0 = 0, ^^^^ 
Sgr0 = -T, Sgrr = "S, Sgge = -r^c'^il. 
I 



From the perturbed Riemann tensor, the perturbed 
Einstein tensor can be expressed in terms of gauge- 
invariant variables by means of 

SGab = SRab - ^RSgAB - ^9abSR, (A6) 
where the perturbed Ricci scalar is 

6R = g^^SR^B + Sg^^RBn. (A7) 
After some (tedious) calculations one gets 



E + n 



:^e-ry^J 3*" + + n" - 2—deT' + 



3 , 



n' ~ 3—a'dgT \ + G^, (* - S) 



(A8) 
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for the purely brane part, while the mixed ones read 

1 



5G, 



5G, 



2^+7+2 



(A9) 



and the purely bulk components are 
1 



1 , i 



SGrr = ^e'^O (3* + 2^dl<^ + 2^a'deT - 3(t' + 2 
2 r r 

,5Gee = ir^c'^c-'^n (3* + S) - 2r2e''^'" + r^c^'a' (S' - 5*') + ^r^e'^ (4a" + Sa'^) (S - f]) 



1 



re 



1 , 7' 



- i (4a" + 5a'2) T. 



(AlO) 



r 



where □ is the brane d'Alembertian defined above and to first order in metric and field perturbations, using 
[Eq. (|85|l ]. Eq. 173|) . this perturbed matter Lagrangian reads 



3. Stress-energy tensor 

In terms of the underlying fields, the stress-energy ten- 
sor stemming from Eq. © reads 

T.s = \ (I^.*)^ (I?B$) + \ {V,<^)^ (2?,$) ^^^^^ 

with £mat given by Eq. to zeroth order, this is given 
by 



c"''' O' 1 c^^ O'^ 

q 2 

2 2 

,A"-fA"t e-T „2'^ + '^^ 



A" - A-t dy A" -I- A-t 



2i 



dip 2 



1 e~^ Q'^ 1 1 c~^ 
'Cmat = --;r^ Y-oV''^-o—V^Q'^-^i'P). (A12) sor is therefore given by 

I 



(A13) 

The purely brane part of the perturbed stress-energy ten- 



e-^ Q' 
J.2 q 



1 e-^ Q'2 1 



2 r2 g2 



1 e-'' g'2 1 e-'' 
2~lf^2~ 



A- - A-t dF A" A-t 



2i dv3 2 



1 e--^ Q'2 1 



2 r2 g2 



1 e-"^ 



2 r 



2 



y(^) 



while the mixed components are 



(A14) 



-1 O' 



A'-f A-t 



j,2 q 



ST^e = aJ ^ (6, - 9') + ^Q^^^ - W'QQ [, 



(A15) 
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and the bulk ones 



X' + X' 



"t 



.x + 



X'^ e""^ 

^<<3(556 



X-X'^ dV X + X'^ 
2i 2 



mi 



STre = 



'le-y 2 
2—^ 

|/2 



1 e 
2~r- 

X' 



,2 ^2 



2 r- 



r 



X + X^ X-. 



— re'- — 

dip 



dVX^X^ 



X'-X'^ ,^X-X^ 



X + x^ 



2 r2 g2 



2 



(A16) 



4. Faraday tensor 

In order to directly derive the perturbed Maxwell equa- 
tions from Eq. l[T!?jl . we have used the following perturbed 
Faraday tensor whose purely brane components vanish 



SHf,, = = <5iJ^", 

and with the mixed parts 

5H„r - (6, - e') , 

-e-^77^''9, (9, - 9') , 
(9e - deO) , 
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^ fir 



(A18) 



(9e - 9e9) . 



The only nonvanishing purely bulk components end up 
being 



6Hr9 = 9'« - doQr 



9L 

q 



(A19) 



Owing to these formulas, one can calculate the perturbed 
Faraday tensor divergence involved in Eq. H19|l by means 
of 



% = dJH^^ + T\JH"^, (A20) 



where we have used the antisymmetry property of SH^b ■ 
The perturbed left-hand side of Eq. (|19|l can be, in turn, 
expressed in terms of the perturbed matter fields by 
means of Eqs. I@J and H73() to give the perturbed Maxwell 
equations ^ to (|5^ . 



APPENDIX B: THE 6D CYLINDRICAL 
SOLUTION 

In this appendix, we closely follow the analysis of 
Ref. [73, generalized to the six-dimensional case, and 
show that the only gravity confining solution to Einstein 



equations in the presence of a negative cosmological con- 
stant is the one used throughout this article. 

We start with the most general static, cylindrically 
symmetric line element for a hyperstring inside which 
one assumes also rotation invariance, namely 

ds^ = g^{r)dt^-dr^-g,{r) {dx^ + dy^ + dz^)^g^{r)de^ . 

(Bl) 

Note that if one also demands Lorentz invariance along 
the vortex world sheet, then one would restrict attention 
to the subset of solutions for which ^^(r) = g^tir), if it 
exists. 

Setting = g'^g^Qs — — det{gAB)i Einstein equations 
take the form 



-gl + Au = 0, 
9i 



(B2) 



where a prime indicates a derivative with respect to r, 
and 



9i92 



9i93 9293 



9i 



4A = 0, (B3) 



which can be combined to yield the simple equation for 
the metric determinant, namely 



u" + -Au = 0. 



With A < 0, the general solution of Eq. (|B4|) is 



u = u+e''/'''^ + u.e"''/'''^ 



(B4) 



(B5) 



with 



-5A/2. We shall come back shortly to this 



solution. 

Eq. (|B4|I can be integrated, with the result that 



(B6) 



where K is an arbitrary real constant, i.e. K"^ > 0. This 
latter requirement stems from the fact that, in order to 
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have actual cylindrical coordinates with < 9 < 2n and 
to avoid a singularity at the symmetry point r — 0, one 
must in general impose that limr-^o u{r) = 0. 

Expanding and integrating Eqs. (|B2|I lead to the solu- 
tion 



(B7) 



where the otherwise arbitrary constants Ki are related 
through 

3K, +K,+K,^0, 
coming from the definition of u in terms of gi, and 



5' 



which is nothing but a rewriting of Eq. (|B3|I . 

Now, with u{0) = 0, the solution (|B5p or a direct inte- 
gration of HB6p gives 



is thus only asymptotically valid. Note that this is in fact 
not problematic in the framework of the Abelian Higgs 
vortex since the point at r = cannot be described by 
the vacuum Einstein equations. 

An interesting point concerning the finite- volume met- 
ric (|B1() is that, seen from far away from the vortex, the 
warp factor can be interpreted as a missing angle, just like 
in the simpler case of a Nambu-Goto string in Minkowski 
space. However, in the case at hand, the requirements of 
both cylindrical and maximal symmetries now demands 
a missing angle of 27r, as can be seen by evaluating the 
diameter of a circle r =const., at a distance r — > oo and 
noting that this diameter vanishes with the warp factor. 

One can also note that changing the radial coordinates 
into 



r = r\ exp 



2rA 



rx 



u — Kr\ sinh 



and therefore the metric functions read 



JO) 



where the gj"'' are three constants of integrations satisfy- 
ing 





2 




sinh ( — ) 




tanh ( I 


VaJ. 







(B8) 



(B9) 
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(0) 



2K^ 
'5A 



and [u'(0)] = K^. Eq. HB9|) is but the 6 dimensional 
generalization for nonvanishing cosmological constant of 
the Kasner metric |75||, |72| already obtained in Ref. [t^ ■ 
For the purpose of confining gravity, one needs a space 
with finite transverse volume = J d9dry/—g = 

2tt Ju(r)dr, i.e. one demands that u{r) goes asymp- 
totically to zero faster than 1/r. The solution ljB8|) 
leads however to = 2TrK cosh{r /ta)]^ , which is un- 
bounded unless K —> 0, in which case it is undefined. 
The general solution is thus useless as one needs to im- 
pose K — 0. However, with = in Eq. (jB5|l . the 
transverse volume remains finite, so gravity ends up be- 
ing actually confined on the hyperstring core. In this 
case, since u(0) ^ 0, the solution is singular at r = 0, 
and is thus incomplete, being unable to describe the vor- 
tex location itself. The warped metric 



de^ 



dr^ 



(BIO) 



with = 2rA = ^10/|A|, transforms the metric into 



- \ ^ - 2 

^ j (Ty^.dx'^dx") - r^de^ - (y) df2 



ds'^ 



(BU) 

With these new coordinates, the hyperstring is located at 
f = fA, and is a compact coordinate, varying between 
and 27r/fA. 

Finally, it is worth pointing out that when a time de- 
pendence is allowed for in either l|B10|l or (|B11|I . and 
metric variables are assumed separable, with gi{r) 
gi(r)Ili{t), then one still has a valid solution provided 
Hj = H3 — C{t — to)"^ = n~^, with C and to arbitrary 
constants of integrations. This solution, which is also 
of constant scalar curvature R — 3A, corresponds, to a 
shrinking extra dimension (a circle located at a given co- 
ordinate distance from the vortex gets smaller with time) 
and an expanding brane interior. As discussed above, 
in the framework of the present article, this solution is 
only asymptotically valid. Nevertheless, according to the 
background field solutions (see Sect. IIV( l. this shrinking 
might be associated with a growing conical singularity in 
the vortex core, i.e. with a time-dependent v^. It is a 
peculiarity of this number of dimensions that the scale 
factor in the brane in this model evolves as ex t^, which, 
recalling the time coordinate to be conformal time, is the 
expected evolution of a matter-dominated universe. 
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